A bilateral shear layer between two parallel Couette flows 
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We consider a shear layer of a kind not previously studied to our knowledge. Contrary to the 
classical free shear layer, the width of the shear zone does not vary in the streamwise direction but 
rather exhibits a lateral variation. Based on some simplifying assumptions, an analytic solution has 
been derived for the new shear layer. These assumptions have been justified by a comparison with 
numerical solutions of the full Navier-Stokes equations, which accord with the analytical solution to 
better than 1% in the entire domain. An explicit formula is found for the width of the shear zone 
as a function of wall-normal coordinate. This width is independent of wall velocities in the laminar 
regime. Preliminary results for a co-current laminar-turbulent shear layer in the same geometry are 
also presented. Shear-layer instabilities were then developed and resulted in an unsteady mixing 
zone at the interface between the two co-current streams. 

PACS numbers: 47.15.St, 47.15.Rq, 47.27.ek, 47.15.-x 



I. INTRODUCTION 

The laminar boundary layer which forms in the shear 
region between two semi-infinite uniform parallel streams 
admits a Blasius type of similarity solution. Lock [l| con- 
sidered streams of two different fluids and with different 
velocities, say U2 and U\, as depicted in figure []Ta). His 
numerical solutions are reproduced in textbooks on vis- 
cous flows, e.g. Panton Q and White In the special 
case where the two fluids are the same Lock [l[ found 
that the solution depended only on the ratio U2/U1 of 
the velocities of the two streams. A particular feature 
of the plane shear layer (also referred to as the classical 
mixing layer) is that the width of the shear zone thickens 
as (xv/Ui) 1 / 2 , where v is the kinematic viscosity and x is 
the streamwise distance from the trailing end of the plate 
where the two semi-infinite streams merge. Klcmp and 
Acrivos U considered the non-uniqueness of this bound- 
ary layer problem due to an indeterminacy of the stream- 
line separating the two streams. Even by including all 
higher-order effects in their analysis, the position of the 
dividing streamline remained indeterminate. 

A special feature of the plane shear layer is that stabil- 
analyses, for instance that by Betchov and Szewczyk 
show that the flow is unstable at all Reynolds num- 
bers. Bhattacharya et al. @ ascribed this peculiar- 
ity to the parallel-flow assumption on which the Orr- 
Sommcrfeld equation is based. They therefore formu- 
lated a non-parallel stability problem and found a criti- 
cal Reynolds number close to 30, below which the flow is 
convectively stable. 

The stability of the flow in a circular shear zone was in- 
vestigated by Rabaud and Couder 0] . They developed an 
experimental apparatus in which the fluid was enclosed 
in a very short and broad cylinder, the top and bottom 
of which were both formed of disks rotating at a certain 
angular velocity surrounded by an annulus rotating at a 
different rate of rotation. In §4.3 of Ref. Q they consid- 
ered a linear model problem in which a planar shear zone 
was formed in the center region of a rectangular duct of 



which the two halves moved in opposite directions. 

In this paper another type of laminar shear layer will 
be considered, namely the plane shear layer which forms 
in the interaction zone between two co-current and fully- 
developed plane Couette flow^j]. Contrary to the classical 
free shear layer [H-dj], the width of the new shear layer 
is independent of the streamwise position. As we shall 
see, however, the extent of the shear zone varies in the 
direction normal to the streams and thus shares some 
features with the linear model problem solved earlier by 
Rabaud and Couder Q . In the present study we further- 
more aim to compare the analytically derived solution 
with a numerical solution of the full three-dimensional 
Navier-Stokes equations in order to justify the inherent 
assumption of unidirectional flow. 



II. MATHEMATICAL FORMULATION 

Let us consider the shear layer which develops be- 
tween two streams with velocities U2 and U\ in the 
(x, y) — plane. While U2 and U\ are taken as constants in 
the classical shear layer theory, let us now assume that 
both U2 and U\ vary linearly with the lateral position. 

We will assume that the flow field is steady and fully 
developed in the streamwise direction, i.e. the three ve- 
locity components U, V , and W are functions only of 
the two coordinates y and z. The Navier-Stokes equa- 



1 The flow beneath a flat-bottomed ship, e.g. a tanker or bulk 
carrier, which travels with a small underkeel clearance to a flat 
sea floor can be described as a turbulent Couette flow at full scale 
but may be laminar at model scale in a towing tankQ|. Two co- 
current Couette flows, as considered here, therefore mimic the 
interaction of the flows beneath two adjacent ships in side-by- 
side operation. 
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respectively. Explicitly, the boundary conditions are 

U(x,y,z = 0) =0, 
U(x,y > 0,z = 2h) =U + , 
U(x,y < 0,z = 2h) =U- 

The velocities far away from the shear zone will therefore 
tend to the linear Couette flow profiles U\ = zU-/2h and 
U-2 = zll + /2h when we assume the flow to be steady and 
fully developed. 
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FIG. 1. (Color online) A shear layer between two parallel 
streams, (a) Generic sketch of a planar shear layer between 
streams with velocities Ui and U2; and (b) two co-current 
planar Couette flows with wall speeds U- and U+. The di- 
mensions of the flow configuration refer to those used in the 
numerical calculations. 



tions therefore simplify to: 



Vd v U + W8 z U = v{d;U + d % z U); (la) 
VdyV + Wd z V = -p-'dyP + v(p 2 yV + d 2 V); (lb) 
p~ x d,P + v(d 2 y W + d 2 W); (lc) 

(Id) 



Vd y W + Wd z W 



d y V 



d z W = 0. 



From the above system of partial differential equations, 
it is observed that the 'secondary motion', i.e. the flow 
(V, W) in the cross-sectional plane, is independent of the 
streamwise velocity component U. This flow is governed 
by eqs. (|lb[) — (lld[) for which the zero solution V = W = 
P = is a valid and consistent solution. If so, equation 
(flaj) for the a;— component of the velocity vector simplifies 
to the Laplace equation: 



d 2 U = 0. 



(2) 



A flow like this can be realized in the interaction zone 
between two co-current plane Couette flows, as shown in 
figure QJb). The lower wall at z = is fixed whereas the 
upper wall at z = 2h moves in the positive x— direction 
with a speed which is discontinuous at y = such that 
the wall velocity is U- and U+ for y < and y > 0, 



III. NUMERICAL SIMULATION 

Let us first solve the full Navier-Stokes equations in 
three-dimensional space and time for a flow configuration 
as shown in figure [ljb). The Reynolds number i?e_ = 
U- h/2v based on half of the speed i/_ of the fastest 
moving wall and half of the wall distance 2h was taken 
as 260 and the wall-speed ratio was U- /U+ = 2.0. This 
Reynolds number is well below the subcritical transition 
Reynolds number below which the plane Couette flow is 
known to be stable. 

The Navier-Stokes equations for an incompressible and 
isothermal flow are solved using a parallel Finite Vol- 
ume code called MGLET 0. The code uses staggered 
Cartesian grid arrangements. Spatial discretization of 
the convective and diffusive fluxes are carried out using 
a 2 nd -order central-differencing scheme. The momentum 
equations are advanced in time by a fractional time step- 
ping using a 2 nd -order explicit Adams-Bashforth scheme. 
The Poisson equation for the pressure is solved by a full 
multi-grid method based on pointwise velocity-pressure 
iterations. The computational grid is divided into an ar- 
bitrary number of subgrids that are treated as dependent 
grid blocks in parallel processing. In the present study, 
the size of the computational domain L x x L y x L z and 
the number of grid points in each coordinate direction 
are 50.24/i x 16. Bh x 2h and 256 x 256 x 64, i.e. com- 
parable to the turbulent plane Couette flow simulation 
of Bech et al. [l0| . Uniform grid spacing is adopted in 
the streamwise and the spanwise directions, while a non- 
uniform mesh is used in the wall-normal direction. Peri- 
odic boundary conditions are employed in the streamwise 
and spanwise directions. No-slip and impermeability con- 
ditions are imposed on the walls. 

The numerical solution of the unsteady Navier-Stokes 
equations converged to a steady state. The present com- 
putations thus yielded a flow field which is steady and 
fully developed in the streamwise direction. This steady 
flow field is illustrated in figure [2] where the stream- 
wise velocity profiles are plotted along the span. It can 
be observed that in regions away from the shear zone 
the streamwise velocity varies monotonically in the wall- 
normal direction in a linear fashion similar to laminar 
plane Couette flows. The shear layer width is minimal 
just below the moving plates where the velocity gradient 
(dU/dy) is high. Due to viscous diffusion the steepness 
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FIG. 2. Spanwise variations of the streamwise velocity U/U- 
at different wall- normal positions z/h. Numerical solution of 
the full Navier-Stokes equation. 
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FIG. 3. Normalized (a) wall-normal velocity W/U- and (b) 
spanwise velocity V/U- contours in the mid-plane. Numerical 
solution of the full Navier-Stokes equation. 



of the velocity profiles is gradually reduced towards the 
stationary wall. In figure [3] non-dimensionalized wall- 
normal velocity W and spanwise velocity V at the chan- 
nel mid-plane are shown. These results clearly illustrate 
that these two velocity components are practically zero 
in magnitude and hence negligible. 



IV. ANALYTICAL SOLUTIONS 

The results from the numerical solution of the full 
Navier-Stokes equation showed beyond any doubt that 
the secondary velocity components V and W are negli- 
gible for all practical purposes. The streamwise veloc- 
ity component U is therefore governed by the Laplace 



equation ([2]) for which an analytical solution now will 
be sought. Assuming separation of variables U(y, z) = 
Y(y)Z(z) we obtain the decoupled equations 



Y" 



-k 2 Y: Z" = k 2 Z 



(3) 



where k is a constant. The general solution to the equa- 
tion for Y in terms of cos(fcy) and sin(fcy) and periodic 
boundary conditions Y(y + 2L) = Y(y) imply a dis- 
cretization of k according to k = nir/L where n S N. 
Standard theory of Fourier series to satisfy all boundary 
conditions gives the answer 



U phc (y,z) 
(U) 



z 

2h 



2A 



E 



sinh 



(2«+l)7r Z • {2n+l)ny 
L bH1 L 



tt — sin h ( 2 "+!) 2 ^ 



2n + 1 



(4) 



where the average of the two wall velocities is (U) = 
t;(U- + U + ). One notes that the velocity, when reduced 
by the average of the two wall velocities, depends on £7_ 
and U+ exclusively through a single parameter 



A=(U+-U-)/(U) 
according to (here and henceforth (U) ^ 0) 
U(y,z) 



(U) 



2h 



Af hjL (y,z). 



(5) 



(6) 



The first term is obviously the simple linear Couette flow 
profile for a single wall moving at velocity (U) . All infor- 
mation about the shear layer is now contained in the 
function fh,L(y,z). We will call the second term the 
'asymmetry term' of U{y, z). 

From Eq. ^ we may obtain the solution for the case 
of Neumann boundary conditions at y = ±oo by taking 
the limit L — » oo. The sum then becomes an integral ac- 
cording to the Eulcr-Maclaurin formula (e.g. Ref 



11 



It is possible to solve the resulting integral in closed 
form, however, the Neumann solution is more elegantly 
obtained using complex analysis. 

The boundary value problem is conformally invariant 
under the Schwarz-Christoffcl-typc map 

4:/i 7T77 

n = y + iz = — artanhC; C = tanh — -. (7) 

The boundaries of the original problem are now mapped 
onto the real ( axis so that z = maps to ( € (—1,1), 
the boundary z = 2h, y > where the velocity is U+ 
is mapped to £ € (1, oo) and z = 2h, y < where the 
velocity is U- to ( e (— 1,— oo). All other points in 
the original strip domain are mapped conformally to the 
upper half of the complex 77 plane. Seeking 10(77) so that 
U{rj) = Im[zo(?7)], we can instead find a function w(C,{ii)) 
satisfying the boundary conditions. Such is the case for 

w(C) = — log(C+l)-% log(l-C) (8) 



4 




FIG. 4. (Color online) Streamwise velocity U(y, z)/U- varia- 
tions along the wall-normal direction z/h. Full numerical so- 
lution (circles) compared with analytical Eq. ((5J) (solid lines) 
at spanwise positions (top to bottom) y/h = —4.2, —2.1, 
-0.95, -0.42, -0.16, -0.03, 0.03, 0.16, 0.42, 0.95, 2.1, 4.2. 
The top two and bottom two lines nearly coincide. 



and with some manipulation we find 

U Nbc (y,z) 



(U) 



z 
2h 



A 

— arctan 

7T 



tanh ( — 
4ft. 



/"KZ\ 

tan — 
\4hJ 



(9) 

It is easy to see that simple linear Couette profiles are 
obtained far aside of the shear layer as they should, not- 
ing that for y ^> h, tanh(^) — > Sg(y), i.e., the signum 
function. 



V. DISCUSSION 

The closed-form analytical solution (|9|) is applicable 
when the flow is subjected to Neumann boundary condi- 
tions in the spanwise direction, whereas the solution (|4|) 
applies for spanwise periodicity We first verified that the 
latter series solution with 2L = I6.8/1 coincided with the 
closed- form solution © in the shear zone. A comparison 
between the latter and the full Navier-Stokes solution is 
provided in figure [4] Recall that the width of the com- 
putational domain on which the Navier-Stokes equations 
were integrated was L y = 16.8ft. whereas the displayed 
velocity profiles only span the central 8.4ft of the shear 
zone. The numerical and analytical solutions agree to 
a relative error of order 0.01%, except in the immediate 
vicinity of the discontinuity in the wall velocity where 
deviations of about 1% are found. Hence the shear layer 
is excellently described by the analytical expression given 
as the second term on the right hand side of Eq. ©, and 
an analytical expression for its width variations in the 
wall normal direction can be found as follows. 

Let us define the shear layer symmetrically about y = 
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FIG. 5. (Color online) Numerical shear layer thickness for 
S — 0.01 (solid line) compared with approximation Eq. (|10p 
(dashed line). 



0, that is, it is the area within which the asymmetry 
term of U (y, z) in Eq. ^ differs from its asymptotic form 
far from the shear region by more than a small relative 
measure 6, typically 1%. Far from the shear layer the 
asymmetry term tends to that for simple Couette flow, 
Sg(y)Az/4ft, so according to equation (T9J) the edges of 
this shear layer are given as the positive and negative y 
solution of 



A 

— arctan 

7T 



/ Try \ (ixz 
tanh — tan — 
\4hJ Uft 



Sg(y)(l-*)^. 



It is clear that in the laminar regime the thickness of the 
shear layer is completely independent of the velocities U— 
and U + . The shear layer is quite obviously symmetric 
around y = such as we have defined it. We take the 
tangent of either side and expand to linear order in 5. 
Noting that for |tanhx| close to unity we may use the 
approximation, tanhx w Sg(x)[l — 2exp(— 2x)}, we find 
the approximate shear layer thickness di ay er as a function 
of z as 



^laycr _ _ 2 / TTZ S 

h ~ 7T ° S \4ft sin §f 



(10) 



For small S this is an excellent approximation except in 
the immediate vicinity of the discontinuity as shown in 
figure 

According to its definition, the shear layer thickness 
vanishes at z = 2ft, but not near the stationary wall at 
2 = 0. This is a consequence of the definition of the 
shear layer, which is the region within which the rela- 
tive difference between the real velocity profile and the 
asymptotic (plane Couette) flow at infinity is above some 
threshold. Thus, although all absolute velocities tend to 
zero at z = 0, the values of U taken at different y rela- 
tive to each other remain finite and nonzero in this limit, 
hence the shape of the shear layer as shown in Fig. [5l 

With our ratio L/2h = 8.4 there is very little mutual 
influence between the shear profile at y — and that at 
y = ±L appearing due to the periodic boundary con- 
ditions. To quantify we plot the difference between the 
periodic and Neumann boundary condition velocity pro- 
files, Eqs. (j4]) and (j9]), respectively, in Fig. HI For compar- 
ison the same is plotted when the period L is halved and 



5 




-4 -2 2 4 

y/2h 



FIG. 6. (Color online) Absolute difference between U in the 
cases of periodic (period 2L) and Neumann [L = oo) bound- 
ary conditions, whose analytical expressions are Eqs ((4]) and 
([9]), respectively. This quantifies the interaction of shear lay- 
ers at y = and y = ±L. We compare the period used in this 
paper, 2L = I6.8/1, with half and twice this period. For the 
longest period the difference is at the level of the cutoff error 
from including just 80 terms in the Fourier series, Eq. Q, 
giving rise to the oscillating behavior in this case 



doubled. For our choice of L/h ratio the influence of the 
periodic boundary conditions remains below the 1% level 
all the way out to the half-way points y = ±L/2. Dou- 
bling the spanwise period further reduces the influence 
by five orders of magnitude. 



level of uj z is reached in the vicinity of the moving wall at 
y/h = 0. In the present flow the variations of ui y and lo z 
are governed solely by viscous diffusion. It is noteworthy 
that in order to maintain u) x = 0, tilting of lo v and uj z by 
means of velocity strains, i.e. uj v dU/dy and u) z dU/dz, 
respectively, are exactly outweighed. 

In order to check the stability of the present flow, the 
numerical integration of the full Navicr-Stokes equation 
was repeated with random noise superimposed on the 
initial flow field. The solution eventually evolved to the 
same steady state as before, thereby suggesting that the 
laminar flow is stable at the Reynolds number Re- = 260 
considered here. 

The stability of the shear-layer at a higher Reynolds 
number was further investigated by Narasimhamurthy et 
al. [IH, where Re- was increased to 1300 while keep- 
ing the i?e+ equal to 260. The computational domain 
and the number of grid points were the same as before. 
This resulted in a velocity ratio of U-/U+ = 5.0. The 
higher i?e_ considered in [13| is well above the subcritical 
transition Reynolds number of 300-3 70 ( see Schneider et 
al. and Tuckerman & Barkley !15]) for which the 

plane Couette flow is fully turbulent. Thereby, the in- 
terface between a laminar and a turbulent plane Couette 
flow was studied in Narasimhamurthy et al. [l3[ rather 
than the present shear layer between two laminar Cou- 
ette flows. The higher Re-, i.e., turbulence, had a dis- 
tinct effect on the interface where shear-layer instabili- 
ties were developed resulting in large-scale interactions 
between the turbulent and the nominally non-turbulent 
part flow. Such a large-scale mixing-zone is shown in fig- 
ure [7] where instantaneous velocity components from 1 1 31 ] 



VI. CONCLUDING REMARKS 

We have considered the shear layer between two co- 
current streams with constant vorticities ai = U+ /2h 
and a\ = U- /2h, respectively. The numerical solution 
of the full 3D Navier-Stokes equation first showed that 
V and W were totally negligible away from the plane of 
the velocity discontinuity. The streamwise momentum 
equation therefore simplified to a Poisson equation ([2]) 
for U(y,z). By assuming spanwise periodicity, the an- 
alytical solution (j4j was derived. As the period tends 
to infinity, the analytic solution simplified to the closed- 
form solution ([9]). 

The width of the shear zone increased monotonically 
from the moving split wall towards the fixed bottom 
plane where the width 2di aycr ks 6.8ft-. Outside of the 
shear region, the two co-current streams are constant- 
vorticity Couette flows where the spanwise vorticity uj y = 
dU/dz is generated by the wall motion and diffused 
downwards. In the shear zone, tu y of the fastest-moving 
fluid is reduced in the lower part and increased as the 
moving wall is approached, whereas u y of the low-speed 
flow is enhanced in the lower part but decreases near the 
upper wall. In addition, the shear zone gives rise to a 
wall-normal vorticity u z — —dU/dy > and the highest 
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FIG. 7. Instantaneous velocity components from the DNS 
data of Narasimhamurthy et al. pl| where U- /U+ = 5.0 and 
Re- — 1300 and Re+ = 260. All profiles are extracted from 
the mid-plane, (a) streamwise velocity U/U- ; (b) wall-normal 
velocity W/U-; (c) spanwise velocity V/U-. 
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are plotted. The secondary-flow in figure [7|c) together 
with the meandering-motion or the unsteadiness in the 
other two velocity components as depicted in figure [TJa, 
b) clearly indicate that the shear-layer is no longer stable 
under turbulent conditions. Thereby, a mixing-layer was 
established in [l3| rather than the present stable shear- 
layer (see [l3j for more details on the momentum trans- 
fer and turbulent diffusion mechanisms occurring in the 
mixing-zone) . 

The dynamics of a laminar-turbulent interface has also 
been investigated in a plane Couette flow configuration 
at Reynolds numbers in the range from 180 to 650 by 
Duguet et al. Contrary to the co-current laminar- 

turbulent Couette flow [la ], initial perturbations were 
introduced into a conventional plane Couette flow driven 
by continuous wall motion. Sufficiently strong perturba- 
tions at sufficiently high Reynolds numbers led to turbu- 
lence localized in the spanwise direction which enabled 
detailed explorations of the interface dynamics. 

The Taylor-Couette flow between rotating cylinders 
approaches the plane Couette flow in the limit of large 
radii and small gaps, as shown for instance by Faisst and 
Eckhardt [l7|. The present flow configuration can thus 
be considered as a limiting case of a Taylor-Couette flow 
with different rotational speeds of the upper and lower 
parts of the driving cylinder. 

Another analogy to the currently considered geome- 
try is that considered experimentally, e.g., by Burin and 
co-workers, in which the end-walls of a Taylor-Couette 
set-up consist of two independently moving rings rotat- 
ing with the inner and outer cylinders, respectively [l8j]. 
Again the present flow configuration is obtained as a lim- 
iting case. 

The shear layer which forms between two co-current 
plane Couette flows shares some similarities with the 
shear zone formed in the junction between the two halves 



of a fluid-filled duct which move in opposite directions. 
Rabaud and Couder Q derived an analytic solution 
for the latter problem by assuming that the secondary 
motion was negligible and thereafter solving Laplace's 
equation ([2]). The width of the resulting shear zone 
was largest midway between the two parallel duct walls, 
whereas the width of the present shear layer increased all 
the way from the splittcd moving plate to the stationary 
wall. Another distinguishing feature of the present flow 
is that linear Couette flow profiles are recovered outside 
of the shear zone such that both spanwise vorticity and 
viscous shear still exist. Outside of the shear zone in the 
flow analysed by Rabaud and Couder Q, the fluid was 
conveyed as a solid body along with the moving halves of 
the duct and neither vorticity nor viscous shear stresses 
prevailed. 

A new shear layer has been introduced in this paper. 
Contrary to the classical shear layer, the width of the 
shear zone varies in the direction perpendicular to the 
shear but is independent of the streamwise direction. 
An analytical solution has been derived which compared 
perfectly well with accurate numerical solutions of the 
three-dimensional Navier-Stokes equations. The solution 
turned out to be independent of the fluid viscosity, which 
implies that the solution is valid for all Reynolds numbers 
sufficiently low for the flow to remain stable. 
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